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ABSTRACT 

R144 is a WN6h star in the 30 Doradus region. It is suspected to be a binary because 
of its high luminosity and its strong X-ray flux, but no periodicity could be established 
so far. Here, we present new X-shooter multi-epoch spectroscopy of R144 obtained at 
the ESO Very Large Telescope (VLT). We detect variability in position and/or shape 
of all the spectral lines. We measure radial velocity variations with an amplitude 
larger than 250 kms -1 in Niv and Nv lines. Furthermore, the Nm and Nv line 
Doppler shifts are anti-correlated and the Niv lines show a double-peaked profile on 
six of our seven epochs. We thus conclude that R144 is a double-lined spectroscopic 
binary. Possible orbital periods range from 2 to 6 months, although a period up to 
one year is allowed if the orbit is highly eccentric. We estimate the spectral types of 
the components to be WN5-6h and WN6-7h, respectively. The high luminosity of the 
system (logLboi/L© ~ 6.8 ) suggests a present-day total mass content in the range of 
about 200 to 300 M , depending on the evolutionary stage of the components. This 
makes R144 the most massive binary identified so far, with a total mass content at 
birth possibly as large as 400 M . We briefly discuss the presence of such a massive 
object 60 pc away from the R136 cluster core in the context of star formation and 
stellar dynamics. 

Key words: stars: individual (RMC 144) - stars: early-type - stars: Wolf-Rayet - 
binaries: close - binaries: spectroscopic - stars: formation 



1 INTRODUCTION 

Massive stars strongly contribute to the light from distant 
star-forming regions and dominate the production of ioniz- 
ing radiation. The most massive stars are important contrib- 
utors to the feedback from stars on the interstellar medium 
(Crowther 20071. Yet, how massive a star can be is still 



poorly constrained. In the quest for a firm grip on the max- 
imum stellar mass limit, hydrogen-rich Wolf-Rayet (WNh) 
stars have been identified as the most massive stellar ob- 
jects ( de Koter et al.|1998 l, with typical dynamical masses 
over 75 M Q ( |Rauw et al.|2004| |Gamen et al.|2006| |Nieme la 



Based on observations collected at the European Southern Ob- 
servatory (ESO) under GTO program IDs 085.D-0704 and 086.D- 



0446. 
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ct al. 2008; Sch nurr et al. 2008 2009). |c7o"wther et al.| ( |20T0 l 
have recently re-analysed known WNh stars in R136, the 
central cluster in 30 Dor, and derived probable present-day 
masses up to 250 M©. Other very massive WNh stars are 
found in apparent isolation in the 30 Dor region, among 
them R145 flSchnurr et al.||2"009| ) and VFTS 682 ( jBesten 



lehner et al.||2011|. Unless these stars are all ejected from 



the cluster core through dynamical interactions or supernova 
explosions, their presence challenges the view in which mas- 
sive stars formed along with a cortege of lower mass stars in 
a dense cluster environment. 

With v = 11.15 | |Breysacher et al.|1999[ ), the WN6h star 
R144 (RMC 144, Brey 89, BAT99-118, HD 38282) in 30 Dor 
is the visually brightest Wolf-Rayet star in the Large Mag- 
ellanic Cloud (LMC) and it is one of the early-type objects 
isolated from the R136 core. Because of its brightness and 
its strong X-ray emission, R144 has been proposed as a bi- 
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Figure 1. Optical spectrum of R144 at the epoch of largest RV separation (#2 in TablefTjl. The main spectral lines are identified. 



nary candidate, but no definitive proof of the orbital motion 
has been obtained so far. Moffat ( 1989 1 performed the first 



multi-epoch spectroscopic campaign on R144 and, in the ab- 
sence of variability, concluded that the object was likely sin- 
Using a more extended observational data set, (Schnurr 



et al. ( 2008 2009 1 detected variability in the He II14686 line 



as well as in broadband polarimetry, although no periodic- 
ity could be established. In this work, we report on the de- 
tection of R144 as a double-lined spectroscopic binary and 
we suggest that the two components are WNh stars with 
a present-day total mass of up to 300 Mq for the binary 
system. 



2 OBSERVATIONS AND DATA REDUCTION 

Intermediate-resolution spectra of R144 were collected from 
April 2010 to April 2011 using the X-shooter spectrograph 



( Vernetetal.|2011| > mounted on VLT/UT2 at ESO's Paranal 
observatory. In total R144 was observed at seven epochs, al- 
lowing us to explore time scales of days, weeks and months. 
X-shooter provides us with complete optical and near- 
infrared (NIR) wavelength coverage from the U to the K- 
band by combining three instrument arms, each operating 
as a separate spectrograph. Given the brightness of R144 
and our scientific aims, the narrowest slit of each arm was 
used, providing a resolving power of 9100, 17400 and 11000 
in the UVB, VIS and NIR arms, respectively. The observa- 
tions were carried out in nodding mode with a nod throw of 
5" along the parallactic angle. Exposures of 90s (resp. 60s) 
for the UVB and VIS arms (resp. NIR arm) were obtained 
at each nodding position. Depending on the weather condi- 
tions, two to four nodding cycles were performed, yielding a 
typical signal-to-noise ratio per resolution element (S/N) of 
over 200 in the UVB and VIS arms and over 100 in the NIR 
arm. The journal of the observations is given in Table [l] 

The data were reduced and flux calibrated using the 
X-shooter pipeline running under the Reflex workflow envi- 



Table 1. Journal of the observations. Cols. 3 to 5 give the 
primary (vi) and secondary (^2) radial velocities measured 
from Nvl4944 (vi only) and Nivl4058 (l>i and v 2 ). 

# HJD a - Nvl4944 Nivl4058 

2,400,000 vi (kms" 1 ) vi (kms" 1 ) v 2 (kms^ 1 ) 



1 


55,309.478 


357 


336 


56 


2 


55,452.888 


119 


380 


65 


3 


55,580.589 


151 


64 


319 


1 


55,585.650 


137 


57 


319 


5 


55,589.681 


159 


61 


308 


6 


55,604.525 


387 


360 


59 


7 


55,673.479 


360 


342 


73 



NOTES - a. Heliocentric Julian date at the beginning of 
the exposure. 



ronment ( Modigliani et aL]|2010 ). The spectra of each arm 
were then normalised to the continuum by fitting a polyno- 
mial to continuum regions. An example of the obtained spec- 
tra is presented in Fig. [l] The R144 spectrum is dominated 
by strong emission lines of hydrogen, helium and nitrogen 
that are characteristic of WNh stars (Sm ith et al.|[l996| |. 
Hel lines show clear P-Cygni profiles. We also detect Nv 
P-Cygni lines at 4604 and 4620 A as well as a single-lined 
N V14944 emission. 



3 VARIABILITY AND RV MEASUREMENTS 

We search for line profile variability in the time series formed 
by our spectra by visual inspection and by using a tempo- 



ral variance spectral analysis (TVS, Fullerton et al. 19961 



Significant line profile variability is detected in most of the 
spectral lines. The associated TVS spectrum is often double- 
peaked, which is indicative of large radial velocity (RV) vari- 
ations, hence of a potential binary motion. We will return 
to this aspect later. 
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Figure 2. Nvl4944, Nivl4058 and Nml4640-42 lines at different 
epochs. The vertical dashed lines indicate the rest wavelength of 
N V, Niv and of both Nm components, shifted to the LMC rest- 
frame (i.e., shifted by +270 kms -1 ). Red lines show our best-fit 
RV measurements from Table [T] 



Fig. [2] shows a selection of Nm, Niv and N v lines ob- 
served during our campaign. Nvl4944 has a single compo- 
nent whose position varies by over 4 A.N IV14058 has single 
comparable components from each of the two stars, hence 
the line is double-peaked with the two peaks changing in 
position. Nml4640-42 effectively has two components from 
each star, thus a total of four. The line profile changes from 
a single peak to double as a result of line blending. 

As illustrated in Fig. [2] the N v line and the dominant 
Nm components present anti-correlated RV variations. Fur- 
ther visual inspection of our spectra indicates that the po- 
sitions of the Si IV and Nv lines also vary in opposite di- 
rection. Finally, the widths and centroids of hydrogen and 
helium lines are found to be variable, but the lines remain 
single-peaked, probably due to their intrinsically broad na- 
ture. 

To quantify the amplitude of the RV variations, we use 
the fitting method described in |Sana et" al, (2013), which ad- 
justs Gaussian(s) to line profiles accounting simultaneously 
for all the observations. We separately adjust the Nvl4944 
and N IV14058 line profiles and we weight the fit by the S/N. 
The best fit Gaussian(s) are overlaid in Fig. [2] The associ- 
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Figure 3. N V14944 and Nivl4058 RV measurements. The dotted 
line indicates the average velocity of stars in 30 Dor. 



ated RV measurements are listed in Table [T] and displayed in 
Fig. [3] A local renormalisation has been applied for N V14944 
given that the line is positioned in the wing of Hell4922. 

The Nvl4944 line is well fitted by a single Gaussian 
and displays a peak-to-peak RV variation of 280 kms . 
The typical measurement errors are about 6 kms -1 . In the 
following, we arbitrarily refer to the component associated 
with the N V14944 emission as the primary, hotter star in the 
system. 

The wings of the Nivl4058 line are not properly rep- 
resented by a Gaussian profile and we limit the fit to the 
central part of the line only. Two Gaussians are used and 
Fig. [2] only shows the sum of the fitted Gaussians. The 
N IV14058 RVs somewhat depend on the wavelength range 
adopted for the fit. Combined with the fact that this line 
profile is not Gaussian, we estimate error bars of the order 
of 15-20 km s -1 on the RV values. While the peak of the blue 
component of N IV14058 remains more or less at the same lo- 
cation, its width and the shape of its blue wing changes. The 
best fitting solution indicates thus that the two components 
of N IV14058 do actually exchange position at epochs #3 to 
5 with respect to the other epochs. In this case, the pri- 
mary and secondary components display a peak-to-peak RV 
variation of about 320 and 260 kms -1 , respectively. 

Overall there is a reasonable agreement between the 
primary RVs obtained from the two lines. Niv!4058 seems 
to display a slightly larger RV amplitude than Nvl4944 al- 
though the obtained values remain within the uncertainties 
of the RV measurement method. 



4 THE NATURE OF R144 
4.1 Spectral type 

We first classify the composite spectrum of R144 as if it were 



a single star. The criteria of Smith et al. (1996), which are 



based on ratios of peak intensities of H, He, C and N lines, 
unambiguously confirm the WN6h type of the composite 
spectrum. Given the blue-shifted weak absorptions in each of 
the He I and He 11 emission lines, the final composite spectral 
type is WN6ha. 

As mentioned earlier, the dominant component of the 
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Nnil4640-42 lines moves opposite to Nvl4944 (Fig. [2). The 
fainter Niii14510-14-18 multiplet shows a similar behaviour. 
This suggests that the Nm lines are predominantly associ- 
ated with the secondary, cooler star. In order to qualitatively 
reproduce the profile variations of the N nil4640-42 complex, 
the primary, hotter star is still expected to make a small 
contribution to Niii14640-42, possibly with a peak intensity 
no larger than half that of the secondary. A peak intensity 
of a third of that of the secondary, or less, would imply 
a WN5h classification, which we consider to be within un- 
certainties. A complete spectral disentangling with a larger 
data set is, however, needed to properly quantify the amount 
of N ml4640 in the primary star. 

If most of the N ill emission in the R144 spectrum orig- 
inates in the secondary star, accounting for the dilution of 
the secondary Nm line by the primary light suggests that 
the secondary may be of a slightly later spectral type. We 
thus adopt WN5-6h + WN6-7h as our best estimate of the 
system spectral types. 



4.2 Orbital constraints 

Mass ratio 

The mass ratio of a double-lined binary system can in princi- 
ple be constrained directly from the RV measurements of the 
two components. We use the orthogonal least-square regres- 



sion method of Sana et al. ( 2006 1 to fit the v\ vs. V2 relation, 



allowing for different apparent systemic velocities for the two 
components. We obtain a mass ratio of M2/M1 = 1.17±0.06, 
which implies that the primary, hotter star is the less mas- 
sive star in the system. While it is possible that systematics 
in the RV measurements of the N IV14058 line might actually 
invert the mass ratio, our results clearly indicate that the 
masses of the two components do not differ by more than a 
few tens of percent from one another. 



Orbital period 

Our data only cover one epoch of large RV variations and 
do not allow us to constrain the periodicity of the system. 
The large RV change (> 200 kms -1 ) between observations 
5 and 6, separated by only two weeks, suggests either a rel- 
atively short period or a significant eccentricity, possibly 
both. Combining the epoch of variability observed in our 



X-shooter data with those in Schnurr et al. (2008 2009) 



does allow us to obtain a first constraint on the orbital pe- 
riod. Significant variability is observed by [Schnurr e t al.| at 
t « 2, 452, 630 and at ti ~ 2, 453, 000, hence a At of about 
370 d. Our RV variations lie at t « 2,455,586 (Table [TJ, 
hence almost exactly 8 x At from epochs to (resp. 7 x At 
from ti). This strongly suggests that the maximum orbital 
period of R144 is about 370 d. 

Submultiple integer values (i.e., P ~ 370/n with n = 
1,2,3,...) constitute a family of possible solutions, where 
n = 4, 5 (P ~ 90 or 70 d) and n > 7 (P < 50 d) are unlikely, 
as large RV variations are predicted at epochs where Moffat 
(19891, Schnurr et al. and/or our X-shooter observations 



10 to 20 days. The most favoured periods are thus about 2, 
4, 6 and 12 months. 

The semi-amplitudes of the observed RV variations (Ta- 
ble III provide a lower limit on the true semi-amplitudes of 
the RV curves. For given stellar masses and eccentricity, they 
constrain the maximum orbital period. Adopting masses of 
120 and 140 M Q for both components (see Sect. |5~lj ), and 
eccentricities of 0.3, 0.5 and 0.7, the upper limit on the pe- 
riod is 4.5, 6 and 10.5 months, which matches the range of 
possible periods discussed above. 



5 DISCUSSION 

5.1 Stellar mass content 

Mass estimates can be obtained from the spectral types. 
WN5h stars have typical (present-day) masses in the range 
130 to 250 M , while W N6h and WN7h stars h ave masses 

between 



and 130 M Q (Crowther et al. 2010 and refer- 



show a constant signal. This deduction assumes that the 
peak of the RV variations is concentrated over a time span of 
15 days but the same constraints are obtained by assuming 



ences therein). This suggests that R144 is formed by two 
> 80 M Q stars, making it one of the most massive binary 
systems known. 

The total mass content of R144 can also be estimated 
from its bolometric luminosity (Lboi)- We use the narrow- 
band photometry of Breysacher et al. ( 1999 v — 11.15, b — 
v = —0.11). We estimate the visual extinction in two differ- 
ent ways: (i) we deredden the X-shooter flux-calibrated spec- 
trum to the Rayleigh- Jeans slope, yielding A' v — 1.1±0.2; (ii) 
we use the color excess E(b—v) ~ 0.21 and R' v = 4.12 (corre- 
sponding to Ry — 3.1, Turner|l982 l, resulting in A' v ~ 0.9. 
Both extinction estimates agree within the uncertainties. To 
avoid overestimating the bolometric luminosity, we conser- 
vatively adopt the latter value and obtain an absolute visual 
magnitude M v = —8.2 for the system. The visual magnitude 
of R144 is thus about 1 mag brighter than other WN6h bi- 
naries such as R145 ( M v ~ -7.2; |Schnurr et al.||2009| and 
WR20a (M v = -7.04; |Rauw et al.||2005[ ). Given that both 
R145 and WR20a have a total mass content of 160 Mq or 
larger, the brightness of R144 suggests an even larger mass 
content. 

While the bolometric correction remains uncertain 
without a detailed atmospheric modelling, it is likely in 
the range 4.0 (a typical early O star) to 4.7 mag (see e.g. 
Bestenlehner et al.|2011 |, yielding logLboi/i-© ~ 6.8 to 7.0. 
The total luminosity of R144 is thus very similar to that of 
the other WNh stars in the core of R136 and brighter than 
VFTS682 in the surrounding nebula. 

Adopting logl/boi/£o = 6.8 and a mass-luminosity re- 
lation appropriate for very massive stars (Kohler et al., in 
prep.), we estimate masses of 170 + 205 Mq if the stars are 
on the zero-age main sequence (MS) and of 80 + 95 Mq if 
the stars are close to the MS turn-off, i.e. the TAMS. 



5.2 "Isolated" formation or dynamical ejection 

R144 is located at a projected distance of about 60 pc from 
the central cluster R136. The presence of very massive stars 
outside cluster cores is not predicted by most massive star 
formation theories, which expect that massive stars form 
as part of a cluster ( jZinnecker fc Yorke||2007[ ). R144 may 
thus have been ejected from the cluster core, in which case 
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R144 would have needed a (projected) runaway velocity of 
~ 60 kms~ to travel from R136 to its current location 
in w 10 6 years. Because the two components of R144 are 
both very massive stars (i.e., with a short evolutionary time 
scale), the supernova kick scenario is unlikely. In the run- 
away scenario, R144 then needs to have been ejected from 
R136 through dynamical interactions. 

In dense stellar environments such as R136, the most 
massive stars tend to sink quickly to the cluster center where 
they dynamically interact with each other. During such in- 
teractions, the less massive object is typically ejected (Heg- 



gie 



1975 1 and the standard theory thus requires an even 



more massive object to have caused the ejection of R144. 
Such an object is not seen. In principle, this object might 
already have terminated its evolution despite the young age 
of the region. However, stars of hundreds of solar masses 
are all expected to have a similar life-time of about 2.3 Myr 
(Kohler et al., in prep.). As other very massive stars are 
still present in the core of R136, this scenario would require 
significant fine tuning. 

An alternative scenario for ejecting R144 from R136 
that does not require such an extremely massive object is 
offered by recent theoretical work. |Fujii fc Porteg ics Zwart 
( 2011 1 calculated that the gravothermal collapse of a cluster 



core produces, in the core, a binary that is formed by the 
most massive stars in the cluster. This binary then domi- 
nates the dynamical interaction in the core, frequently eject- 
ing other stars. It hardens by each dynamical encounter until 
it gets ejected as well. In that scenario, R144 is a candi- 
date "bully binary", possibly formed in the gravothermal 
collapse of the cluster core. This hypothesis is not without 
drawbacks. First, single stars more massive than the R144 
components likely exist in the cluster core ( jCrowther et al.| 
20101, so that the bully binary is not formed by the most 



massive stars. Second, its orbital period (< 370 d) seems too 
short compared to the typical periods of systems predicted 
to form during the collapse (P ~ 1000 d, Fujii & Portegies 
Zwart|201l| >. 



If R144 is not a runaway system, it has thus formed in 
situ and in relative isolation from R136. One possibility is 
that R144 is part of an older association along the line of 
sight toward the northern part of 30 Dor. R144 is indeed 
located in the vicinity of the luminous WN stars R146 and 
R147 (projected distance of about 22pc), and within a gi- 
ant high-velocity/X-ray shell that is most likely a supernova 
remnant ( jFeast et aL]|1960| |Walborn et al.||2013| ). This sce- 
nario, however, supposes the entire dispersion of the parent 
association to explain the low surface density of O- and B- 
type stars in the vicinity of R144. 

The last scenario, in which R144 has actually formed in 
situ and independently of a parent cluster or OB association, 
would pose a serious challenge to massive star formation 
scenarios that rely on a dense stellar environment, such as 
formation through mergers of lower-mass proto-stellar cores 



( Bonnell & Bate 2002 1 or through competitive accretion 
( Bonnell fc Bate|2006 l. The formation of very massive stars 
in relative isolation is compatible with a fractal vision of star 
formation, where a giant molecular cloud fragments in vari- 
ous cores which then form stars stochastically. As discussed 
in |Bressert et alT] ( |2012[ ), discriminating between these sce- 
narios is not only important for our understanding of high- 



mass star formation but also to discern how nature samples 
the initial mass function. 



5.3 Evolutionary stage 

While awaiting an accurate determination of the orbital 
properties, the fact that the primary, hotter star seems to 
be less massive than the secondary, cooler component is in- 
triguing. A lower-mass primary component is typically the 
indication of past or ongoing binary interaction. Given the 
relatively large separation of the R144 system, the stars are 
only expected to fill their Roche lobe towards or after the 
end of the MS. It is, however, possible that a high eccentric- 
ity allows for interaction close to periastron passage. The bi- 
nary orbit may also have been widened by mass loss through 
stellar winds or as a result of a past mass-transfer event, so 
that R144 may have been a tighter binary in the past. 

Alternatively, the primary star may have been the more 
massive component at birth. Given its larger mass, it may 
have entered the WR-wind mass-loss regime earlier than the 
secondary, thus shedding mass at such an accelerated rate 
that it became the less massive component in the system. 
A qualitative comparison with the Kohler et al. (in prep.) 
evolutionary tracks indicate that a 260+175 Mq pair at zero- 
age MS evolves, after 2 Myr, into a 90+120 Mq system with 
a total present-day luminosity of 10 6 ' 9 Lq and, indeed, with 
a hotter and less massive primary component. 

Characterizing both the orbit and the physical prop- 
erties of the R144 components may further help to eluci- 
date the evolutionary status of the object, hence to better 
constrain stellar evolution at very high masses. Gathering 
high quality spectra with a sufficiently good phase cover- 
age to enable the determination of the full orbital solution 
and to perform spectral disentangling is thus a high priority 
goal in gaining a better understanding of R144. A linear po- 
larimetry monitoring campaign and a search for atmospheric 
eclipses may further help to constrain the inclination of the 
system (e.g., Moffat et al. 1998), thus to provide a direct 



measurement of the masses of the two components of R144. 
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